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A3STRAC' 


As  a  sequel  to  the  1955  paper  by  Eorrin-Eroukovsfcy  this  study  re¬ 
presents  an  additional  step  in  the  gradual  development  of  the  theory  of 
ship  notions  in  regular  waves  which  was  started  by  Eriioff  in  l£9-  and 
which  has  gained  aoaentus  in  recent  years  in  the  work  of  aeixrlua,  St, 
Denis,  norvin-Xroukovsky  and  Le»~L£.  Certain  errors  and  omissions  of  the 
previous  work  are  corrected  and  the  theory  ir.  its  present  fora  is  veri¬ 
fied  by  comparison  of  the  theoretically  computed  notions  of  eight  widely 
different  ship  models  with  the  notions  observed  in  towlng-tank  tests 
over  a  range  of  node!  speeds  and  at  several  wave  lengths.  An  excellent 
correlation  is  demonstrated  in  the  case  of  a  destroyer  model,  the  shape 
of  which  conforms  most  closely  to  the  assumptions  of  the  theory.  Coed 
agreement  is  found  in  the  cases  of  several  typical  commercial  ship  forms. 
The  agreement  is  generally  satisfactory  'Doth  in  amplitudes  of  heaving  and 
pitching  notion  and  in  the  phase  relationships  between  the  notions  and 

‘The  material  in  this  paper  was  prepared  at  the  Experimental  Towing  Tank, 
Steven s  Institute  of  Technology  as  a  part  of  the  activity  under  Office  of 
Naval  research  Contract  No.  hour  263-10  (ETT  Project  PI  169"?).  repro¬ 
duction  in  whole  or  in  part  is  permitted  for  any  purpose  of  the  !.  S. 
love  rumen  t. 

2 

research  Professor  of  Fluid  Dynamics,  Experimental  Towing  Tank,  Stevens 
Institute  of  Technology,  Hoboken,  N.  J. 

■^Mathematician,  Exoerinental  Towing  Tank,  Stevens  Institute  of  Technology. 

Paper  to  be  presented  at  the  Annual  Meeting,  New  fork,  N.  I ,  November 
195",  of  the  Society  of  N'aval  Architects  and  Marine  Engineers. 


the  waves.  The  theory  is  found  to  fail  only  for  a  sailing  yacht  cha¬ 
racterized  by  large  slopes  of  sides  at  the  UrfL  throughout  its  length 
and  by  a  pronounced  bov  overhang. 

A  discussion  is  given  of  the  usefulness  of  the  theoretically 
computed  responses  to  regular  waves  in  predicting  ship  motions  lx  an 
irregular  sea  aiyi  a  satisfactory  agreement  is  shown  for  the  destroyer 
model  between  the  statistical  averages  of  heaving  and  pitcldng  anpli*’ 
tudes  of  motion  obtained  in  this  way  and  those  obtained  experimentally* 

INTRODUCTION 

In  a  papor  written  for  the  November  1955  meeting  of  the  Society 
of  Naval  Architects  and  Marine  Engineers,  the  first  author  presented 
a  theory  of  the  heaving  and  pitching  motions  of  a  ship  in  regular  head 
or  following  seas*  That  paper  added  one  more  chapter  to  the  historic 
development  which  was  started  by  Kriloff  (I896,  1898)\  and  which 
received  a  strong  impetus  in  recent  years  from  the  work  of  Weinblum 
and  Su*  Denis  (1950),  St*  Denis  (1951)  and  Korvin-Kroukovsky  and 
Lewis  (1955)*  The  particular  advances  made  by  the  first  author  in  his 
1955  paper  were  consideration  of  heave-pitch  coupling  and  a  theoretical 

^Tbe  Bibliography  at  the  end  of  the  paper  is  arranged  alphabetically 
by  authors'  names,  with  the  date  of  publication  given  in  parenthesis,, 
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calculation  of  the  exciting  forces  exerted  on  a  ship  by  wares  the 
results  of  which  were  in  satisfactory  agreement  with  experimental 
measurements.  The  calculation  took  into  account  the  interaction  bet¬ 
ween  ship  and  wave  which  is  neglected  in  the  heretofore  universally 
applied  "Froude-Kriloff  hypothesis". 

The  motions  of  two  ship  models  computed  on  the  basis  of  the  t  h?ory 
were  compared  with  the  motions  measured  in  a  towing  tank.  A  reasonably 
good  agreement  in  the  amplitudes  of  pitching  and  heaving  was  demonstra¬ 
ted,  but  there  was  a  large  systematic  error  in  the  phase  relationships. 

In  the  1955  paper  a  "strip"  method  of  analysis  was  used,  not  only 
because  of  its  mathematical  simplicity,  but  also  because  by  this  method 
the  distribution  of  hydrodynamic  loads  along  the  hull  could  be  obtained. 
This  distribution  combined  with  the  distribution  of  inertial  loads  re¬ 
sulting  from  the  ship  motions  can  be  directly  applied  to  tbs  analytical 
calculation  of  bending  moments  truing  on  a  ship  structure  in  waves. 

Since  publication  of  that  paper,  the  second  author  has  been  engaged  in 
evaluating  ship  bending  moments  on  this  basis"*  for  comparison  with  ex¬ 
perimental  results.  It  was  found  in  the  process  that  bending  moment 

'’Under  the  auspices  of  Panel  J>-3  of  the  Hull  Structure  Conmittee  and 
Panel  H-7  of  the  Hydrodynamics  Committee  of  the  Society  of  Naval 
Architects  and  Marine  Engineers. 
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calculations  require  a  much  greater  precision  than  calculations  of 
motions,  and  that  the  theoretical  method  which  gave  reasonably  good 
amplitudes  of  motion  was  not  yet  adequate  for  the  calculation  of  bend¬ 
ing  moments. 

In  the  course  of  re-examination  of  the  original  method  occasioned 
by  this  difficulty,  several  errors  and  omissions  were  discovered  and 
rectified.  Most  important  has  been  the  realization  that  the  velocity- 
dependent  (or  damping)  terms  of  the  equations  of  motion  consist  not 
only  of  energy  dissipative  parts  but  of  dynamic  non-dissipative  parts 
as  well.  The  existence  of  the  latter  had  been  demonstrated  by  Haskind 
(19U6)  and  Havelock  (1955),  but  not  sufficient  attention  had  been  paid 
to  their  studies,  perhaps  because  of  the  difficulty  in  following  the 
mathematics  of  the  first  and  because  of  the  relatively  recent  appear¬ 
ance  of  the  second. 

Introduction  of  the  necessary  corrections  has  improved  the  corre¬ 
lation  between  the  calculated  and  experimental  amplitudes  of  motion 
and  has  resulted  in  a  generally  good  correlation  of  the  phase  relation¬ 
ships,  Furthermore,  confidence  in  the  method  has  increased  with  more 
extended  application  from  the  two  models  studied  in  1955  to  a  total  of 
eight  models  of  widely  different  types  at  present.  The  pitching  and 
heaving  motions  of  these  models  in  regular  waves  have  been  experimentally 
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measured  and  analytically  calculated,  and  the  results  are  described 
in  this  paper. 

In  his  paper  of  1955  the  first  author  emphasized  that  his  objective 
was  to  place  in  the  hands  of  naval  architects  a  working  tool  to  use  in 
the  process  of  developing  better  ship  forms.  It  is  the  aim  of  the  present 
paper  to  present  an  improved  tool  and  a  broader  basis  for  appraisal  of 
its  success  in  application. 

THE  PATTERN  OF  SEAKEEPING  RESEARCH 

Before  proceeding  with  a  detailed  analysis  of  ship  motions  in  re¬ 
gular  head  seas,  it  may  be  advisable  to  demonstrate  how  this  particular 
subject  fits  into  the  general  pattern  of  research  into  the  seakeeping 
qualities  of  ships.  The  complete  picture  would  require  an  investigation 
of  ship  motions  in  all  six  modes  —  surge,  sway,  heave,  roll,  pitch  and 
yaw  —  in  realistic  complex  seas  of  any  direction.  Such  a  complete  in¬ 
vestigation  has  not  yet  been  made,  but  as  a  first  step  the  simpler  pro¬ 
blem  of  ship  motions  in  head  or  following  seas  has  been  tackled  by- 
several  writers,  beginning  with  Kriloff  (1896).  A  further  simplifica¬ 
tion  has  been  to  assume  regular  seas,  so  that  side  sway,  rolling^  and 

^  The  authors  are  indebted  to  Prof.  Edward  V.  Lewis,  Mr.  Edward  Numata 
and  Mr.  John  Dalzell  for  all  experimental  data. 

7 

Grim  (1952)  shows  that  at  certain  speeds  in  a  following  sea  rolling 
instability  develops  despite  all  symmetry. 


■  t 


yawing  notions  are  eliminated  and  only  surging,  heaving  and  pitching 
notions  remain.  Of  these,  sieging  appears  to  have  little  effect  on 
heaving  and  pitching,  so  that  it  is  sufficient  to  consider  the  ructions 
in  these  two  modes.  Only  in  certain  aspects  of  towing-tank  testing  tech¬ 
nique  has  surging  appeared  to  be  important  (Vedeler,  1955  >  Sibul,  1956, 
and  Reiss,  1956). 

Formerly  the  response  of  a  ship  to  regular  wave 3  was  thought  of 
as  an  approximation  to  the  actual  behavior  at  sea.  Since  the  publica¬ 
tion  of  papers  by  St.  Denis  and  Pierson  (1953)  and  by  Fuchs  and  McCamy 
(1953) >  and  the  subsequent  work  of  Korvin-Kroukovsky  and  Lewis  (1955), 
Lewis  (1955),  Lewis  and  Numata  (1957),  Cartwright  and  Rydill  (1957)  and 
Cartwright  (1957),  the  point  of  view  has  changed.  Obtaining  the  res¬ 
ponse  of  a  chip  to  regular  waves  is  no  longer  considered  as  an  end  in 
itself,  but  rather  as  the  end  result  of  the  hydro-mechanical  phase  of 
the  broader  problem,  which  is  to  be  followed  by  a  stati c  •  .  al  phase  in 
order  to  obtain  the  response  to  realistic  complex  (or  irregular)  seas. 

It  should  be  emphasized  that  the  ship  properties  such  as  shape  and  mass 
distribution  are  considered  only  in  the  hydro-mechanical  phase;  the 
statistical  phase  is  a  strictly  mathematical  treatment  of  tas  results 
obtained  in  the  first.  The  results  of  the  first  phase  are  required 
only  to  be  in  the  form  of  a  plot  of  amplitudes  of  pitch  and  heave  per 


unit  wave  height  versus  frequency  of  encounter.  Various  methods  of 
non-diaensional  plotting  proposed  and  used  fros  time  to  tine  hare  little 
usefulness  in  this  connection. 

The  ship  response  to  regular  waves  can  be  measured  in  a  tewing 
tank  vithou.  uuy  reference  to  theory.  However,  experimental  data  do 
not  bring  out  readily  the  significance  of  'the  various  physical  proper¬ 
ties  of  a  ship.  It  is  necessary  to  develop  a  reasonably  comprehensive 
theory  in  order  to  gain  an  understanding  of  the  effects  of  these  pro¬ 
perties  on  the  notions  of  a  ship.  Such  a  theory  would  also  provide  a. 
basis  for  judging  the  experimental  data.  There  are  nary  pitfalls  in 
towing  tank  testing  of  ship  models  in  waves  which  can  be  avoided  when 
knowing  what  to  expect  on  theoretical  grounds. 

Tbe  theoretical  approach  to  the  proolets  of  ship  responses  to 
regular  waves  has  taken  two  general  directions.  One  is  to  undertake 
tbe  complete  problea  of  a  tody  oscillating  in  waves  at  once,  i.  e.  to 
include  initially  in  the  equations  of  notion  the  dynamics  of  the  body 
aad  the  hydrodynamics  of  a  fluid  with  a  free  surface.  Probably  the 
most  complete  study  of  this  kind  is  that  of  Has kind  (19ii6),  which  has 
been  aade  accessible  to  English-speaking  readers  by  Dr.  rfehausen's 
translation  (under  the  auspices  of  the  5—5  Panel  of  tbe  Society  of 
Naval  Architects  and  Marine  Engineers).  These  studies  have  been 


valuable  in  providing  guidance  to  the  simpler  approaches*  but  their 
mathematical  difficulty  and  the  need  for  extensive  calculations  to 
bring  out  the  final  results  prevent  their  being  used  directly  by  naval 
architects*  The  other*  simpler*  approach  has  been  taken  by  Kriloff 
(1896*  1898),  Weihblum  and  St.  Denis  (195'0)  and  St.  Denis  (1951).  Here 
the  differential  equations  of  motion  are  written  on  the  basis  of  the 
dynamics  of  rigid  bodies  and  carious  coefficients  are  independently 
and  individually  evaluated  by  hydrodynamic  considerations.  They  are 
visually  derived  by  analogy  with  a  fully  submerged  body  and  the  effects 
of  a  free  water  surface  are  allowed  for  by  separate  considerations. 
Specialized  investigations  of  virtual  masses  and  damping  such  as  those 
by  Holstein  (1936,  1937  a*b),  Havelock  (I9l2),  Ursell  (19ii9,  1953*  195U), 
and  Grim  (1953)  are  available,  and  while  these  studies  are  by  no  means 
uninvolved  the  results  are  presented  in  simple  form  and  can  be  readily 
included  in  the  solution  of  the  equations  of  motion.  With  this  sup¬ 
porting  literature  the  simpler  theoretical  approach  has  become  more  and 
more  realistic. 

The  present  paper  follows  this  simpler  approach.  Its  immediate 
objective  is  the  calculation  of  ship  responses  to  regular  head  or  fol¬ 
lowing  seas*  with  the  expectation  that  this  material  can  then  be  statis¬ 
tically  treated  to  permit  prediction  of  ship  behavior  in  a  realistic 
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irregular  sea.  In  addition  it  is  desired  to  provide  a  method  that  can 
be  readily  extended  to  the  calculation  of  ship  bending  moments. 

THE  EQUATIONS  OF  HOTION 


Neglecting  the  surging*  the  equations  representing  the  equilibrium 
of  forces  and  moments  acting  on  a  ship  moving  in  regular  head  or  follow¬ 
ing  seas  can  be  written  in  the  form 


m  z  +  c  s  "  Hj, 

J  9  «■  C  9  *  H 

m 


(1) 


where  m  is  the  mass  of  the  ship  itself  and  J  is  its  moment  of 


inertia;  z  is  the  heaving  displacement  and  9  the  angular  displace¬ 
ment  in  pitch;  c  and  C  are  the  unit  restoring  force  and  moment 
respectively  due  to  ch.-vvrr-s  in  the  buoyant  forces  arising  from  the 
deviation  of  a  ship  from  its  attitude  of  normal  flotation;  is 
the  force  and  the  moment  about  the  center  of  gravity  of  a  ship  due 

to  water  pressures  generated  in  the  flow  by  the  action  of  the  wave  as 


well  as  by  the  ship's  own  oscillations.  Ttese  forces  and  moments  are 


computed  in  the  Appendix*  and  are  found  to  be  expressible  by  polyno¬ 
mials  the  terms  of  which  are  divided  into  those  corresponding  to  the 
ship's  oscillation  in  smooth  water  and  those  corresponding  to  direct 
wave  action  on  a  restrained  ship.  When  the  terms  of  the  first  kind  are 


transferred  to  the  L.H.  side  of  Equations  (l),  the  equations  take  the 


form 

a*z+bz  +  cz  +  d(?+e8  +  g9a5r  ^ 
lIS  +  BO  +  Ce  +  Ds+Es+Gz  =•  M  eiu)t 


(2) 


The  complete  expressions  for  the  coefficient  on  the  L.H,  side  are 
given  in  the  Appendix  by  Equations  ( U2 )  •  The  terms  on  the  R.H.  sides 
of  Equations  (2)  are  known  as  the  exciting  force  and  moment.  They  are 
given  in  complex  form  in  order  to  facilitate  the  algebraic  work  of  the 
solution  and  it  Is  understood  that  only  the  real  part  of  the  exponential 
on  the  R.H.  side  of  Equations  (2)  is  to  be  taken.  In  this  form 

F  *  F^e*0"  the  real  part  of  F  »  Fq  cos  (a)  t  +  0-), 

_  i  i  /i\+  *  (3) 

M  *  M  e  and  the  real  part  of  H  e  =*  M  cos  (uj  t 
o  o 

Here  Fq  and  are  respectively  the  amplitudes  of  heaving  force  and 

pitching  moment  caused  by  waves,  as  computed  by  Equations  (26)  and  (27) 

in  the  Appendix;  and  X  are  the  phase  angles  at  which  the  maxima  of  the 

force  and  moment  occur  (positive  angles  are  leads,  negative  angles  are 

lags).  (T  or  T  is  zero  if  the  maximum  occurs  at  the  instant  when  the  wave 

crest  is  amidships.  The  harmonic  oscillation  of  force  and  moment  is 

given  by  the  factor  exp(ipit)  in  Equations  (2),  where  cO,  written  for 

brevity  in.  place  of  u)  ,  is  the  frequency  of  wave  encounter.  The  terms 

in  9  in  the  first  and  in  z  in  the  second  of  Equations  (2)  represent 


the  interaction  between  hearing  and  pitching  motions.  These  terms  are 
usually  referred  to  as  "cross-coupling  terms". 


It  is  assumed  that  the  ship  is  moving  in  a  uniform  sea,  and  that 
a  steady  state  of  heaving  and  pitching  is  established.  In  such  a  case, 
the  transient  responses  have  been  damped  out,  and  only  a  particular 
solution  of  Equations  (2)  is  required.  Since  the  exciting  functions 
are  assumed  to  be  harmonic,  the  solution  is  likewise  harmonic,  i.  e. 
of  the  form 

z  =  Z  e'^*’  and  0  =  0  e^1^  (U) 

Z  and  0  are  complex  amplitudes 

Z  =  Z  e1<S  and  0  -  0  e16 

o  o 

where  is  the  absolute  value  of  the  heaving  amplitude,  0Q  that  of 

pitching,  and  -  <5  'cid  -  £  are  the  phase  lags  of  ship  motion  relative 
to  t. 


The  solution  of  Equations  (2)  is  readily  obtained  as 

w  HQ-FS 
4  'QE-PS 


-  FR-KP 
b  '  Q  R  -  P  S 


(5) 

(6) 


where  the  capital  letters  represent  the  groupings  of  the  coefficients 


of  Equations  (2)  as  follows* 
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P  **  -a  +  i  b  ui  ♦  c 

p 

5*-du)  ■fieu)'(’g 

L  (7. 

R«-Du)~*iEu)"-G 

S  *  -Au)2  M  Bu\  +  C 

were  the  " cross-coupling*  coefficients  of  Equations  (2)  set  equal  to 
zero,  i.e.  d“e*g*D*E“G"0>  the  two  equations  would  represent 
two  simple  harmonic  oscillators  jin  forced  motion,  and  Equations  (5)  and 
(6)  would  take  the  familiar  form  of  "magnification  factor11. 

THE  SHIP  MODELS  USED  FOE  COMPUTATIONS  AND  EXPERIMENTS 

The  heaving  and  pitching  motions  in  regular  head  seas  have  been 
rale -.Hated  from  Equations  (5),  (6)  and  (7)  above  for  eight  ship  models 
of  different  forms.  The  properties  of  these  models  are  given  in  Table  1 
and  the  body  plans  in  Figs.  1  to  5.  These  hulls  are  currently  subjects 
of  research  at  the  Experimental  Tewing  Tajik  of  Stevens  Institute  of 
Technology  for  various  projects  sponsored  severally  by  the  David  Taylor 
Model  Basin,  the  Office  of  Naval  Research,  and  the  Society  of  Naval 
Architects  and  Marine  Engineer's.  This  collection  of  models  does  not 
therefore  represent  a  series  developed  for  a  particular  purpose;  it 
includes  a  wide  range  of  ship  forms  and  displacement-length  ratios-. 

These  forms  may,  however,  be  considered  in  groups  of  two  or  three  having 
certain  resemblances  and  differences. 
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Hodel  IhhS  is  the  Series  60,  0.60  block  coefficient  fora  and  Model 
l6l6  is  an  experimental  modification  of  it  by  £.  ?.  Levis  (1955)  which 
has  the  sane  dine  ns  ions  and  the  sane  afterbody  but  an  extreae  7-type 
forebody  (Fig.  l).  These  two  models  were  investigated  by  Korvin-Kroukovsky 
in  his  1955  paper  but  have  been  submitted  to  re-analysis  by  the  corrected 
and  improved  procedure  given  in  the  present  paper.  The  newly-computed 
values  of  all  the  coefficients  for  wave  length  X  equal  to  ship  length 
L,  X/L  *  1,  are  listed  in  Tables  2  and  3,  which  replace  Tables  h  and  6 
of  the  earlier  paper,  and  the  new  curves  of  B,  b,  and  versus 
frequency  and  of  C  and  g  versus  forward  speed  are  shown  in  Figs.  6 
and  7  which  replace  Figs.  1  to  U  in  the  earlier  paper.  The  calculated 
and  experimentally  measured  resnnTtsea  of  these  rod  els  to  waves  of  two 

lengths,  X/L  =  1  and  1.5,  are  shown  in  Fig.  8  which  replaces  Figs.  5 
to  8  in  the  earlier  paper. 

Model  lliiiii  is  the  T-2  Tanker  (Fig.  2)  which  has  been  the  subject 
of  many  experiment-?  conducted  by  E.  V.  Lewis  (195k)  in  connection  with 
the  measurement  of  bending  moments  (under  the  sponsorship  of  the  S-3 
Panel  of  the  Society  of  Naval  Architects  and  Marine  Engineers).  The 
calculated  and  experimen tally  measured  notions  for  X/L  “  1  are  shown 


in  Fig.  9, 
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Model  1723  is  the  destroyer  PD692  (Lewis  and  Dalsell,  1957)  of 
displacement-length  ratio  60  (Pig.  3).  The  results  of  computations  and 
towing  tank  tests  for  four  different  wave  lengths  are  shown  in  Pig.  10. 

Model  1699A  (Fig.  U)  is  the  German  fishing  trawler  "Stralsund* , 
described  by  Captain  W.  MMekel  (1953)  as  possessing  particularly  good 
seagoing  qualities.  It  is  a  short  stubby  hull  of  displ aceaent-iength 
ratio  2 Oh.  Model  I699C  has  the  same  section  forms  as  the  trawler  but 
spaced  wider  apart  and  dimensioned  to  give  the  low  displacement-length 
ratio  of  60  like  the  destroyer.  Pig.  11  shows  the  calculated  and  ex¬ 
perimentally  measured  motions  of  1699 A  and  1699C  in  waves  of  two  lengths, 

X/L  —  1  and  1.25.  Attention  should  be  called  to  the  different  speed 
scales  on  the  L.K.  and  R.H.  sides  of  Pig.  11. 

Model  1699B  (Fig.  5)  is  the  British  cruising  yacht  "Brsmbling" , 
also  reputed  to  have  excellent  sea-going  characteristics  (Fox,  1936). 

With  a  displacement-length  ratio  of  370,  it  is  the  stubbiest  model  of 
the  group  investigated.  The  lengthened  counterpart  of  the  yacht. 

Model  1699D,  has  the  same  section  forms  spaced  wider  apart  and  a  dis¬ 
placement-length  ratio  of  60  like  the  destroyer  and  the  lengthened 
trsxler.  Pig.  12  compares  the  calculated  and  experimentally  measured 

1  ^  f 

1 


actions  of  1695®  and  1699®  in  waves  of  two  lengths,  X/L  ”  1  arac 
X/L  *  1.25.  (liote  difference  in  speed  scales  on  L.H.  and  2.H.  sices 
of  this  figareO 

Since  stach  snail  vessels  as  trawlers  sad  yachts  have  succeeded  in 
operating  efficiently  in  an  open  ocean,  it  was  thought  that  their  feres 
aay  have  seme  particularly  desirable  features  worth  investigating. 

E.  7.  Lewis  (1955)  and  Senate  and  Lewis  (1957)  have  -v  vr.  the  advantages 
of  a  low  displ scener t-le ngth  ratio  so  that  the  high  displacement-length, 
ratio  of  these  saaT'l  ships  appears  to  be  a  drawback.  Tee  reason  for  the 
seaworthiness  of  these  balls  nest  lie,  then,  in  their  freeboards  and 
oectiny,al  for*?. 

The  destroyer  (Model  1723),  the  lengthened  trawler  (Model  1695*0) 
and  the  lengthened  ya'.nt  (Model  1699®)  can  be  compared  for  the  effect 
of  section  fora  at  the  sane  low  cispiaceaert-lesgtfa  ratio  of  fO.  The 
effect  of  section  fora  is  also  shown  an  ?ig.  3  by  oaaparison  of  Model 
ILii 5 ,  the  Series  60  hull,  and  Model  1616.  Model  1616  has  the  saae 
sectional  area  cmrre  and  the  saw  afterbody  as  Model  ILLS  btrt  its  foie— 
body  is  notified  to  hare  extaeae  7  sections,  acre  or  less  1  Ve  a  Maier 
fora. 

Tie  effect  of  greatly  reducing  the  displacesaent-length  r  ...o. 
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1.  e.  of  stretching  the  models,  while  keeping  the  same  section  foras 
can  be  seen  by  comparison  of  the  original  and  the  lengthened  versions 
of  the  trawler  and  the  yacht  on  Figs.  11  and  12. 

CORRELATION  OF  CALCULATED  AND  EXPERIMENTAL  RESULTS 

Excellent  agreement  between  analytically  calculated  and  experiment¬ 
ally  measured  results  is  found  on  Fig.  10  which  depicts  the  responses  of 
the  destroyer  (Model  1723)  to  waves  of  four  lengths,  X/L  -  1,  1.25,  1.5 
and  2,  which  cover  the  most  inport  ant  operational  range  o  It  should  be 
noted  that  of  the  eight  this  model  comes  closest  to  being  "wall-sided" 
at  the  LWL,  and  in  this  sense  conforms  closely  to  the  assumptions  made 
in  the  theoretical  derivations  of  the  Appendix  as  well  as  in  the  support¬ 
ing  literature  on  damping. 

Next  in  degree  of  correlation  co*.i  me  Series  60  hull  (Model  liili5). 
Fig.  8,  and  the  T-2  tanker  (Model  lUUi),  Fig.  9.  Here,  too,  the  agree¬ 
ment  is  excellent  between  the  calculated  and  experimentally  measured 
pitching  motions  and  ^hase  relationships.  However,  there  is  a  signifi¬ 
cant  overes tination  of  the  heaving  motions  at  and  near  synchronism,  for 
both  models,  and,  while  the  calculated  phase  lags  in  heave  are  close  to 
the  measured  lags  for  Model  lUt5,  for  Model  lJqlih  there  is  soiae  discre¬ 
pancy  at  speeds  higher  than  2  ft/sec.  Both  of  these  models  have  U  sections 
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forward,  where  they  deviate  little  from  the  assumption  of  vall-sidedness. 
Pronounced  slopes  of  the  tangents  to  the  section  form"  at  the  UrfL  occur 
only  at  the  stem  where  the  relative  ship-water  motions  are  generally 
smaller  and  the  effect  of  section  form  is  less  important. 

The  calculated  and  observed  amplitudes  of  motion  for  Model  1616 
(lower  half  of  Fig.  8)  agree  very  well  at  X/L  =  1,  but  at  X/L  =1.5 
there  is  a  significant  error  in  the  estimate  of  pitching  amplitude.  The 
phase  relationships  show  good  correlation  in  pitch  but  not  in  heave. 

This  model,  with  its  extreme  V  sections  at  the  bow,  has,  of  course,  a 
pronounced  slope  of  the  tangents  to  the  sections  at  the  1ML,  deviating 
in  this  respect  from  the  assumptions  of  the  theory. 

In  the  case  of  the  lengthened  trawler  Model  1699C  (right  hand  of 
Fig.  ll),  there  is  on  the  average  a  very  good  agreement  between  the 
calculated  and  test  results,  but  the  calculations  underestimate  the 
pitching  anplitudes  and  do  not  show  the  exaggerated  narrow  peak  in 
heaving  at  synchronism.  The  calculated  phase  relationships  agree  very 
well  with  the  test  results.  As  can  be  seen  from  Fig.  h ,  the  slopes  of 
the  tangents  to  the  bow  sections  are  more  moderate  for  the  trawler  than 
for  the  u-bov  Model  1616  but  greater  than  for  Kcael  lUhS- 

The  original  trawler  Model  16 99A  combines  these  section  character¬ 
istics  with  the  high  displacement-length  ratio  of  2Ch.  The  experimental 
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and  calculated  notions  agree  quite  well  at  X/L  ■  1,  but  there  are  large 
discrepancies  in  the  pitching  amplitudes  and  phase  lags  at  X/L  *  1.25. 
The  agreement  in  heave  amplitudes  and  phases  at  the  latter  wave  length 
is  excellent. 

The  only  models  for  which  the  computation  method  fails  completely 
are  the  yachts,  both  the  original  Model  1699B  and  the  lengthened  Model 
IA99D.  These  hulls  have  not  only  large  slopes  of  the  tangents  to  the 
sections  at  the  UfL  at  bow  and  stern,  but  also  appreciable  slope  at 
midihip  iiction  and,  in  addition,  a  largi  bow  overhang  and  a  cutaway 
forefoot.  These  deviaticus  from  the  assumptions  of  the  present  theory 
of  ship  motions  appear  sufficient  to  mutilate  computed  results. 

It  v&*  expected  that  estimates  of  the  motions  of  the  original 
yacht  would  be  poor  not  only  because  of  the  slopes  of  the  sides  at  the 
UfL  and  the  large  overhang  but  also  because  of  its  s tubbiness.  This 
model  was  included  in  the  analysis  in  an  attempt  to  establish  the  limits 
of  the  applicability  of  the  theory.  Ifewever,  the  failure  in  the  case  of 
the  lengthened  yacht  with  a  low  displacement-length  ratio  of  60  is  illu¬ 
minating.  It  indicates  that  the  section  slopes  and  the  end  overhang 
have  much  greater  weight  in  the  calculations  than  tbs  displacement- 
length  ratio.  The  xelatively  small  importance  of  the  latter  could  have 
been  inferred  frca  Vosser’s  (1956)  calculations  of  three-dimensional 


dancing  forces  for  three  length-beam  ratios,  as  shown  in  Fig.  17,  and 
from  the  comparison  by  Korvin-Kroulcovsky  (1955a)  of  his  inertial  force 
calculations  for  a  spheroid  of  length-diameter  ratio  5  with  Havelock's 
for  an  infinitely  long  spheroid. 

To  summarize,  there  is  excellent  agreement  between  the  calculated 
and  experimentally  measured  ship  motions  in  the  case  of  a  destroyer, 
and  a  generally  satisfactory  agreement  in  the  case  of  normal  ship  forms 
Only  in  the  case  of  a  yacht  is  the  computation  method  completely  in¬ 
adequate. 

In  the  discussion  above  the  authors  have  purposely  refrained  from 
mentioning  non-line arity.  Certainly  the  slopes  of  the  hull  surface  at 
the  ML  indicate  a  strong  non-linearity  in  the  restoring  forces  and  the 
corresponding  cross-coupling  terns,  i.  e,  in  the  coefficients  c,  C,  g 
and  G  of  the  equations  of  motion.  However  it  is  not  certain  whether 
it  is  this  non-line rrity  or  the  inexactness  of  the  virtual  mass  and 
damping  coefficients  which  is  responsible  for  the  discrepancies  between 
the  calculated  and  the  observed  motions.  Attention  should  be  called  to 
the  fact  that  although  the  fores  for  which  Grim  has  derived  the  damping 
forces  include  what  appear  to  be  extreme  V  sections,  these  sections 
have  a  sharp  turn  of  the  bilge  and  are  tangent  to  a  vertical  at  the 


ML,  whereas  the  bow  sections  of  a  Maier  form  or  of  the  yacht  discussed 


above  have  large  slopes  at  the  JjfL  and  so  may  well  have  quite  different 
damping  properties.  Unfortunately  neither  theoretical  nor  experimental 
information  on  such  sections  is  available.  Information  available  on 
wide  V  sections  used  in  speed  boats  and  in  seaplane  hulls  has  been  derived 
for  very  high  Froude  numbers,  at  which  the  force  of  gravity  can  be  ne¬ 
glected  in  comparison  with  acceleration  forces,  and  is  not  valid  for  the 
range  of  Fronde  numbers  found  in  ship  operation. 

THE  SIGNIFICANCE  OF  THB  CHARACTERISTIC  PROPERTIES  OF  SHIP  FORMS 

As  a  result  of  all  investigations  made  so  far,  the  effects  of 
certain  properties  of  ship  forms  have  become  established.  The  natural 
pitching  or  heaving  period  of  a  ship  oscillating  in  smooth  water  becomes 
somewhat  modified  in  coupled  motion,  but  nevertheless  remains  a  mc>3t 
important  criterion.  The  shorter  the  natural  pe-iod,  the  higher  the 
frequency  of  wave  encounter  at  synchronism,  and  therefore  the  higher 
the  8 hip  i  eed  for  synchronism  in  waves  of  any  given  length.  The  waves 
causing  large  chip  motions  are  equal  to  or  longer  than  ship  length,  and 
E.  7.  Lewis  (1955)  has  shown  that  the  limiting  ship  speed  in  rough  wea¬ 
ther  therefore  incieases  with  ship  length.  This  effect  i3  clearly  shown 
on  Figs.  11  and  12  by  a  comparison  of  the  original  and  lengthened  models 

g 

of  the  trawler  and  the  yacht0. 

Q 

Note  that  the  speed  scales  are  different  for  the  original  and  length¬ 
ened  hulls. 


However,  a  naval  architect  will  arrive  at  a  certain  ship  length 
and  a  certain  displacement-length  ratio  by  nan 7  considerations  other 
than  seakeeping  ability.  Whatever  the  dimensions  chosen,  it  is  desir¬ 
able  to  see  what  improvement  in  seakeeping  can  be  obtained  by  select¬ 
ing  a  suitable  ship  fora.  Evidently  the  most  important  characteristic 
here  is  the  damping.  The  theoretical  work  of  Holstein  (1936,  1937a, 
1937b),  Havelock  ( 19it2 )  and  Grim  (1953)  shows  that  damping  increases 
rapidly  with  decrease  of  the  mean  draft  of  a  ship  section.  Yet  avoid¬ 
ance  of  bow  emergence  and  slamming  dictates  a  large  actual  draft.  Both 
conditions  lead  to  the  requirement  of  a  low  section  coefficient, 

usually  associated  with  a  V  form,  or  at  least  with  a  large  deadrise  of 

9 

the  3 hip  bottom  .  This  is  indeed  a  characteristic  of  the  small  ships 
operating  in  high  seas,  the  yachts,  fishing  trawlers  and  coast-guard 
cutters. 


These  principles  are  well  illustrated  by  the  comparison  on  Fig.  8 
of  the  Series  60  Model  IhhS  with  its  modified  Y-bow  version.  Model  l6l60 
The  larger  beams  of  the  bow  sections  of  Model  l6l6  increase  the  moment 


Q 

'  The  advantage  of  a  low  section  coefficient  in  reducing  the  anplitude 
of  heaving  notion  in  waves  is  vividly  demonstrated  in  the  recent  theore¬ 
tical  work  of  Grim  (1957).  Of  two  sections  having  draft  equal  to  half- 
beam,  the  one  with  section  coefficient  of  0.555  is  shown  to  oscillate 
with  a  magnification  factor  close  to  unity  through  a  wide  range  of  fre¬ 
quencies,  while  the  one  with  section  coefficient  of  1.015  has  a  magni¬ 
fication  factor  of  3  at  a  sharp  synchronous  peak. 
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of  inertia  of  the  vaterplane  and  this  results  in  a  higher  ship  speed 
for  synchronism.  The  smaller  mean  drafts  of  the  V  sections  increase 
the  damping  forces  and  thereby  markedly  decrease  the  amplitudes  of 
pitching  and  hearing. 

The  advantages  of  the  trawler  and  yacht  section  forms  over  those 
of  a  conventional  ship  like  the  Series  60  lie  largely  in  the  increased 
damping  due  tc  smaller  section  coefficients.  There  is  rather  little 
difference  in  the  behavior  of  the  destroyer  (Model  1723)  and  the  length¬ 
ened  trawler  and  yacht  (Models  1699C  and  1699B),  which  were  made  to  match 
the  displacement-length  ratio  of  the  destroyer,  because  the  destroyer 
also  has  tls*  desirable  low  section  coefficients.  Indeed  the  very  low 
meaa  draft  of  the  stern  sections  of  the  destroyer  gives  it  somewhat 
heavier  damping  and  slightly  lower  amplitudes  of  motion. 

This  advantage  of  the  destroyer  over  the  lengthened  trawler  is 
offset  in  part  by  the  fact  noted  during  the  model  tests  (Nurnata  and 
Lewis,  195?)  that  the  destroyer  frequently  shipped  water  more  heavily. 

It  appears  that  concentration  of  heavy  damping  at  the  stern  is  not  as 
desirable  as  a  more  uniform  distribution  of  damping  between  bow  and 
stem.  This  conclusion  is  in  agreement  with  the  often-expressed  opinion 
that  " do uble -ended*  boats  are  preferable  for  good  seakeeping  qualities, 
and  in  fact  trawlers ,  and  fishing  boats  in  general,  show  a  distinct  trend 


in  this  direction* 


K*p™,DOCUmen* 
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The  beneficial  effect  on  the  behavior  of  a  ball  resulting  from  a 
low  section  coefficient  with  its  corollary  a  large  damping  coefficient 
is  quite  evident  qualitatively  from  the  familiar  definition  of  "asspli- 
fication  factor" ,  i.  e.  the  solution  of  the  simple  uncoupled  equations 
of  notion.  It  is  impossible,  however,  to  judge  the  effects  of  the  <rross- 
coupling  coefficients  by  inspection  of  the  coupled  equations  of  notion 
or  their  solution.  Experimentally  also  these  effects  are  difficult  to 
trace  because  of  the  masher  of  different  .odels  required  and  the  tecious- 
ness  of  towing- tank  testing.  So  far  all  analytical  efforts  have  been 
expended  in  developing  a  theory  of  strip  notion  and  in  demons trating  its 
validity.  Until  recently  the  theory  was  not  good  enough  to  warrant  the 
series  of  calculations  necessary  to  bring  out  the  influence  of  the  cross- 
coupling  coefficients;  but.  new  it  appears  that  at  least  for  ship  forms 
not  deviating  too  much  from  conventional  ones  the  theory  is  sufficiently 

*0  »*fc»  BU«h  PftjupuUUeni  ¥<mhwhile,  K  eomputaUnnal  aeries 
of  systematic  changes  of  form  leading  to  changes  of  cross-coupling  coef¬ 
ficients  could  be  easily  organized  and  carried  out  economically  on  many 
of  the  currently  available  computing  nit  .chines. 
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PREDICTION  OF  MOTIONS  IN  IRREGULAR  SEAS 
PROM  THE  CALCULATED  RESPONSES 

As  stated  earlier,  the  ship  responses  to  regular  waves  do  not 
indicate  directly  the  behavior  in  irregular  waves.  However,  a  statis¬ 
tical  treatment  of  these  responses  in  combination  with  the  spectrum  of 
sea  wavfcs  can  give  results  which  are  significant  for  the  ship  motion 
in  the  realistic  sea.  The  technique  of  this  process  was  first  described 
by  St.  Denis  and  Pierson  (1953)  and  it  was  applied  to  model  tests  by 
E.  V.  Lewis  (1955) •  A  simple  exposition  was  presented  by  Korvin-Kroukovsky 
(1956,  1957).  In  this  process  the  square  of  the  amplitude  of  ship  motion 
per  unit  wave  height,  for  a  given  wave  length  at  a  particular  speed  (i.e. 
for  a  given  frequency  of  encounter), is  multiplied  by  the  ordinate  of  the 
sea  spectrum  (corrected  to  that  ship  speed)  at  the  same  frequency  of 
encounter.  These  products  plotted  against  frequencies  of  encounter 
(resulting  from  a  range  of  wave  lengths  at  a  particular  speed)  form  the 
spectrum  of  the  ship  motions  at  that  speed.  Since  the  ship  motions  in 
an  irregular  sea  are  also  irregular,  varying  continuously  in  amplitude 
and  period,  they  can  only  be  described  statistically,  as  for  instance 
by  mean  amplitude,  mean  of  the  1/3  highest  amplitudes,  or  mean  of  the 
1/10  highest  amplitudes.  The  last  figure  is  often  taken  as  the  probable 
near-maximum.  These  means  were  derived  by  I-ongue t -Higgins  (1952}  as 
multiples  of  the  square  root  of  the  area  E  under  the  response  spectrum 
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curve:  the  mean  amplitude  is  equal  to  .88  the  mean  of  the  l/3 
highest  is  equal  to  1.U1  '/_E,  and  the  mean  of  the  l/lO  highest  is  equal 
to  1.8  Hi!.  These  factors  are  applicable  to  the  typical  spectra  of 
ocean  waves  as  well  as  to  the  ship  motions  in  then. 

As  an  illustration  of  the  above,  pitch  and  heave  spectra  were 
predicted  for  the  destroyer  model  No.  1723  from  the  calculated  respon¬ 
ses  to  regular  waves  of  eight  lengths  at  two  ship  speeds.  The  squares 
of  uhe  computed  responses  were  divided  by  the  square  of  the  wave  height 
and  the  quotients  were  multiplied  by  the  ordinates  of  the  irregular  wave 
spectrum  derived  from  analysis  of  the  waves  in  the  towing-tank  tests, 
thus  yielding  motion  spectra.  These  spectra  were  compared  with  the 
spectra  obtained  experimentally  by  Prof,  E.  V.  Lewis  and  his  associates 
at  the  Experimental  Towing  Tank  of  Stevens  Institute  of  Technology 
(Lewis  and  Dalzell,  1957)  who  used  two  different  methods:  first,  the 
spectra  were  obtained  directly  by  statistical  analysis  of  irregulsi- 
model  motions,  and  second,  ths  spectra  were  computed  from  experimentally 
measured  responses  to  a  series  of  regular  waves.  A  comparison  of  the 
vsrioB  apsetra  is  given  in  Fig.  13,  where  the  solid  lines  represent  the 
spectra  resulting  from  the  statistical  analysis  of  ths  experimental 
records  of  ship  motions  in  irregular  waves,  the  dash  lines  the  spectra 


computed  from  the  experimental  responses  to  regular  waves,  and  the 
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crosses  the  spectra  predicted  by  the  methods  of  the  present  paper. 

It  should  be  emphasized  that  the  shape  of  the  spectrum  curve  has 
no  significance  for  the  amplitude  of  ship  motion.  Components  of  all 
frequencies  exist  simultaneously  and  only  the  area  of  the  spectrum,  Z, 
is  significant.  I^eans  of  the  10$  highest  amplitudes,  computed  from 
the  areas  under  the  spectra,  are  given  in  Table  iu  The  amplitudes  of 
heaving  motion  obtained  by  the  three  different  methods  are  in  excel¬ 
lent  agreement.  The  pitch  amplitudes  computed  from  responses  to 
regular  waves  are  slightly  smaller  than  those  resulting  directly  from 
a  statistical  analysis  of  the  irregular  model  motions.  On  the  average 
there  seems  to  be  no  important  difference  between  the  amplitudes  com¬ 
puted  from  the  experimental  and  those  computed  from  the  theoretic?! 
responses  to  regular  waves  (second  and  third  lines  of  Table  U) j  at 
zero  speed  the  theoretically  computed  values  are  nearer  the  statis¬ 
tically  derived  ones,  at  2.53  ft. /sec.  the  values  based  on  experiments 
in  regular  waives  are  closer  by  about  5$.  The  deviations  are  well  with¬ 
in  the  accuracy  of  the  calculation  procedure.  Since  there  is  always  a 
rather  large  statistical  uncertainty  in  evaluating  sea  spectra,  further 
improving  the  accuracy  of  the  calculation  of  responses  to  regular  waves 
would  have  little  meaning.  It  is  only  important  to  look  for  ways  to 
correct  consistent  discrepancies  or  a  trend  contrary  to  reality. 
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The  discussion  so  far  has  been  limited  to  amplitudes  of  motion. 

These  are  readily  obtainable  by  both  experimental  and  theoretical  methods 
in  use  at  present  but  they  are  not  important  just  by  themselves.  Other 
phenomena  may  be  more  decisive  in  qualifying  the  seakindliness  of  ships. 
The  accelerations,  which  for  a  given  frequency  are  proportional  to  the 
amplitudes,  are  more  important  for  passenger  and  transport  ships  since 
there  is  a  direct  connection  between  accelerations  and  sea-sickness 
(Geller,  19U0,  and  Shaw,  195U).  On  fishing  trawlers  the  accelerations 
impose  hardships  on  the  crew  at  work  (Mtfckel,  1953)=  On  modern  cargo 
ships,  with  crew  accommodations  not  far  from  amidships,  the  accelera¬ 
tions  in  pitch  and  hea-we  are  of  minor  Importance  except  for  the  structure 
supporting  the  cargo  in  the  lower  decks  in  No.  1  hold.  The  critical 
conditions  limiting  the  operation  of  these  ships  appear  to  be  shipping 
of  water  and  slamming,  the  latter  mostly  in  the  light  condition.  These 
conditions  probably  also  limit  the  operations  of  such  naval  ships  as 
destroyers. 

The  statistical  treatment  of  the  accelerations  is  similar  to  that 
of  the  amplitudes  of  the  motions.  The  •7ertical  acceleration  at  the 
bow,  for  example,  can  be  darivod  analytically  from  the  computed  motions 
by  the  formula  u>  sq,  where  aQ  is  the  amplitude  of  vertical  acce¬ 

leration,  sq  is  the  amplitude  of  vertical  displacement  obtained  by 
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combining  the  heaving  and  pitching  motions  in  the  correct  phase,  and 
u)  is  the  frequency  of  encounter.  Acceleration  spectra  can  be  obtained 
and  from  these  the  mean  amplitudes,  the  mean  of  the  l/lO  highest  ampli¬ 
tudes,  etc.,  in  irregular  waves. 

E.  V.  Lewis  (1955)  has  shewn  how,  by  superposition  of  the  path  of 
the  ship's  bow  on  the  regular  wave  profile,  the  bow  submergence  or 
emergence  can  be  evaluated  for  each  wave  length  and  ship  speed.  If 
these  values  are  treated  statistic slly  by  the  method  outlined  above 
for  amplitudes  of  motion,  the  frequency  of  deck  submergence  or  fore¬ 
foot  emergence  in  an  irregular  sea  can  be  estimated.  The  frequencies 
thus  obtained  in  one  case  were  found  to  -^ree  well  with  observations 
in  a  towing  tank  irregular  wave  test.  Since  the  ship  actions  derived 
analytically  here  for  hulls  of  normal  fora  are  close  to  the  model  test 
data,  the  calculated  motions  tan  be  used  instead  of  model  tests  to  pre¬ 
dict  bow  submergence,  hence  shipping  of  water,  and  frequency  of  forefoot 
emergence,  which  would  qualitatively  indicate  likelihood  of  slamming. 

It  is  well  known,  however,  that  out  of  a  number  of  forefoot  emer¬ 
gences  only  a  few  result  in  slamming.  In  order  to  predict  quantitatively 
the  frequency  of  slamming  a  more  elaborate  statistical  treatment  of  either 
experimental  or  calculated  responses  to  regular  waves  is  needed.  The 
probability  must  be  found  of  the  joint  occurrence  of  forefoot  emergence, 


a  high  descending  velocity  of  the  bow  and  approximate  tangftncy  of  the 


keel  to  the  water  surface  at  the  instant  of  impact.  A  beginning  has 
been  made  in  the  work  of  L.  J.  Tick  (195U)  where  the  probability  of 
the  joint  occurrence  of  two  of  the  conditions,  forefoot  emergence  and 
high  descending  velocity,  was  investigated. 


It  is  believed  that  the  phase  relationships  are  of  greatest 
importance  in  this  connection.  Sea  captains  and  other  observers  are 
aware  of  the  difference  small  changes  of  speed  make  in  the  shipping  of 
water  and  slanming.  With  propeller  R.P.M.  on  a  typical  cargo  ship  re¬ 
duced  to  say  75  in  mildly  adverse  conditions,  the  change  of  5  R-P.-1- 
often  makes  the  difference  between  acceptable  and  unsatisfactory  ship 


bsnavior.  Tnsrs  is  little  change  in  amplitude  of  motion  in  this 


and  the  change  in  ship  behavio-  appears  to  be  primarily  due  to  change 
in  phase  relationship.  At  the  lower  speed  the  descending  ship  bow  may 
settle  quiotly  into  the  hollow  of  a  wave  and  at  5  R.P.M.  more  it  may 
strike  with  force  the  flank  of  an  oncoming  wave,  causing  either  slam¬ 
ming  or  shipping  of  water  or  both  in  succession. 


Since  the  phase  relationships  in  regular  waves  are  given  Snfcis** 
factorily  by  the  method  of  calculation  presented  hers,  it  can  be 
expected  that  this  behavior  will  be  predicted  correctly.  Thus  the 
stage  has  been  reached  at  which  the  significant  aspects  of  ship  per¬ 
formance  at  sea  can  be  predicted  by  analytical  methods. 
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COW CLUDING  HEK1BKS 


At  the  conclusion  of  lorrin-Kroukorslcy’ s  1955  paper*  two  lines  of 
research  were  suggested  to  improve  the  theory  of  ship  motions  in  regu¬ 
lar  wares*  These  were  l)  the  solution  of  non-linear  equations  of  notion 
(using  numerical  methods  of  computation  if  necessary)  and  2)  the  expe¬ 
rimental  evaluation  of  the  various  coefficients  of  the  equations  of 
motion.  The  need  for  such  research  appeared  to  be  acute  since  there 
was  a  large  systematic  error  in  the  analytically  computed  phase  lags* 
However,  has  been  shown,  most  of  this  error  was  due  to  (mission  of 
certain  dynamic  coupling  terms.  Application  of  the  *-r  -~.ctad  linear 
theory  described  in  the  present  paper  has  yielded  generally  satisfactory 
estimates  of  both  amplitudes  and  phase  relationships  for  ships  of  normal 
commercial  fora.  Apparently  the  deviations  of  most  of  these  ship  foras 
from  the  simplifying  assumptions  of  the  theory  are  not  large  enough  to 
affect  the  results  seriously. 

While  the  suggested  research  is  no  longer  critically  necessary  it 
is  still  desirable.  The  present  theory  has  failed  in  application  to  a 
rather  unusual  ship  fora  (the  yacht  and  its  lengthened  counterpart) 
characterised  by  sloping  sides  at  the  IML  throughout  its  length  and  by 
a  pronounced  bow  overhang  and  a  cutaway  forefoot.  Development  of  the 


theory  along  the  lines  proposed  is  absolutely  necessary  for  success  in 
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estimatlng  the  behavior  of  such  hulls  and  would  make  the  theory  more 
reliable  for  ship  types  like  the  Haier  forn  and  ths  fishing  trawlers. 

Since  publication  of  the  theory  of  ship  notions  in  irregular  seas 
by  St.  Denis  and  Pierson  (1953)*  a  large  aaount  of  practical  experience 
in  its  application  has  been  accumulated  by  E.  V.  Lewis  (1955*  1957)* 
Levis  and  Nunata  (1957)*  and  Lewis  and  D;  Lze 11  (1957)  in  connection 
with  U  wing  tank  experiments,  and  by  Cartwright  and  Rydill  (1957)  and 
Cartwright  (1957)  in  connection  with  an  actual  ship  at  sea.  The  vali¬ 
dity  and  usefulness  of  this  statistical  theory  have  now  been  fully 
established.  The  theory  of  motions  in  regular  waves  can  therefore  no 
longer  be  considered  by  itself  alone.  It  is  only  a  part  of  the  picture* 

o']  nhpieo  which  SStsbllshSS  tbs  Oi*  £  cK^pi  c 

notions  on  its  fora  and  mass  distribution.  The  results  obtained  in 
this  part*  the  ship  responses  to  regular  waves  or  "response  factors", 
are  then  treated  by  methods  of  mathematical  statistics  in  conjunction 
with  a  measured  or  assumed  spectrum  of  a  realistic  irregular  sea  to 
give  the  realistic  ship  motions. 

These  statistical  methods  as suae  that  the  notions  vary  linearly 
with  wave  height  and  require  that  the  results  of  the  hydro-mechanical 
phase  be  supplied  in  linear  fora  as  ratios  of  response  to  wave  height. 
The  rocomnendation  that  solutions  of  non-linear  equations  of  motion  be 
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obtained  should  thus  be  qualified.  What  is  desired  is  not  a  rigorous 
solution  of  a  non-linear  problem  but  rather  a  solution  of  a  substitute 
linear  problem  which  would  approximate  the  true  solution.  Such  a  pro¬ 
cess  of  "equivalent  linearization"  has  been  used  by  Golovato  (1956)  to 
take  into  account  the  effects  of  non-linear  damping  in  heave.  However, 
non-linearities  are  usually  pronounced  in  both  the  damping  force  coef¬ 
ficients  and  the  restoring  force  coefficients,  and  the  coupling  of 
herring  end  pitching  notions  further  aggravates  the  difficulties  of 
this  probl  em. 

It  is  to  be  noted,  also,  that  in  the  statistical  method  the  ship 
responses  to  a  number  of  wave  lengths  are  averaged,  so  that  non- 
eyetemmtio  errors  tend  to  compensate  each  other.  Because  of  this,  and 
the  statistical  uncertainties  in  the  evaluation  of  a  sea  spectrum  as 
well,  there  is  less  need  for  a  high  degree  of  accuracy  in  estimating 
the  ship  response.  The  theory  of  ship  motions  can  thus  in  its  present 
state  of  development  be  practically  useful.  The  application  of  the 
theory  is  not  limited  to  amplitudes  of  motions  and  accelerations;  the 
frequency  of  shipping  water  can  be  readily  predicted  and  prediction  of 


the  frequency  of  slamming  is  not  too  far  in  the  future. 
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NOMENCLATURE 

The  following  nomenclature  is  used  throughout  the  paper* 

*  Coefficients  of  miscellaneous  terms  of  the  differential 
equations  of  motions.  Equations  (2). 

*  Coefficients  of  miscellaneous  terms  of  the  differential 
equations  of  motions,  Equations  (2). 

“  Ratio  of  amplitude  of  waves  made  by  ship  to  amplitude 
of  heaving  motion. 

=  Amplitude  of  vertical  acceleration  at  the  bow. 

*  Beam  (local). 

-  Instantaneous  distances  of  ship  bow  and  stern  from  the 
nodal  point  of  wave  as  defined  in  Fig.  Vi, 

■  Wave  celerity. 

-  Area  under  spectrum. 

*  Hydrodynamic  heaving  force. 


A,  B,  C,) 
D,  E,  G  ) 

,  b,  c,  ) 

>  b  t 

I 

a 

o 

B 

b,  s 

c 

E 

F 


Heaving  force  imposed  on  a  ship  by  waves  (  =  ¥ ^e11^). 
Acceleration  of  gravity. 

Force  due  to  water  pressures  generated  by  waves  and 
the  ship's  oscillations. 

Moment  about  c.g.  due  to  water  pressures  generated  by 
waves  and  ship's  oscillations. 

Wave  amplitude. 

Longitudinal  moment  of  inertia  of  a  ship  in  mass  units. 
Coefficients  of  Equation  (25) . 

Added  mass  coefficient  in  two-dimensional  vertical  flow 
about  a  ship  section. 

Correction  coefficient  for  effect  of  free  water  surface 
Ship  length. 

Hydrodynamic  moment. 

Pitching  moment  imposed  on  a  ship  by  waves  (»  M^e  ). 
Mass  of  a  ship  or  of  a  ship  section. 

Vertical  damping  force  per  unit  of  body  length  per  foot 
per  second. 

Groupings  of  coefficients  of  the  differential  equations 
of  motions  defined  by  Equations  (7). 


Pressure. 


Radial  distance  to  a  point  q  in  fluid. 

Local  radius  of  semi -cylindrical  body. 

Sectional  area. 

Amplitude  of  vertical  displacement  at  the  bow. 

Time. 

Horizontal  component  of  orbital  velocity  of-water  in  waves 
Ship  speed. 

Vertical  velocity. 

Vertical  component  of  wave  orbital  velocity. 

Longitudinal  co-ordinate  with  respect  to  wave  nodal  point. 

Vertical  co-ordinate  or  local  half-breadth  of  LWL  plane. 

Complex  amplitude  of  heaving  motion  (=  Z  ). 

o 

Vertical  co-ordinate  or  heaving  displacement. 

Polar  co-ordinate. 

Angle  between  longitudinal  tangent  to  body  surface  r>nd 
x-axis. 

Phase  angles  of  ship  motions. 

Local  wave-height  co-ordinate. 
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Angle  of  pitch. 


Complex  amplitude  of  pitching  motion  (■  Q^e1^  )• 


■  Ware  length. 


Longitudinal  co-ordinate  with  respect  to  CG. 


*  Water  density. 


<r3r 


Phase  angles  of  exciting  forces. 


Velocity  potential 


tO  or  u)e  3  Frequency  of  wave  encounter. 


APPENDIX 

EVALUATION  OF  HIDRODINIKIC  FORCES 

This  Appendix  is  arranged  to  follow  as  closely  as  possible  Appendix  I 
of  the  earlier  paper  (Korvin-Kroukovsky,  195-5b)  so  that  the  changes  made 
can  be  easily  seen.  Where  sufficient  discussion  was  given  in  that  refer¬ 
ence,  the  detail  ill  be  omitted  here. 

Formulation  of  the  Problem 

Consider  a  ship  moving  with  a  constant  forward  velocity  V  (i.e. 
neglecting  surging  motion)  with  a  train  of  regular  waves  of  celerity  c. 
Assume  the  set  of  coordinate  axes  fixed  in  the  undisturbed  water  surface, 


with  the  origin  instantaneously  located  at  the  wave  nodal  point  prece¬ 
ding  the  wave  rise,  as  shown  in  Fig.  12*  •  With  increase  in  time  t  the 
axes  remain  fixed  in  space,  so  that  the  water  surface  rises  and  falls 
in  relation  to  them.  This  vertical  displacement  at  any  instant  and  at 
any  distance  x  is  designated  ^  .  Imagine  two  control  planes  spaced 
dx  apart  at  a  distance  x  from  the  origin,  and  assume  that  the  ship 
and  water  with  orbital  velocities  of  wave  motion  penetrate  these  control 
surfaces.  Assvrne  that  the  perturbation  velocities  due  to  the  presence 
of  the  body  are  confined  to  the  two-dimensional  flow  between  control 
planes,  i.e.  neglect  the  fore-and-aft  components  of  the  perturbation 
velocities  due  to  the  body,  as  in  the  "slender  body  theory"  of  aero¬ 
dynamics.  This  form  of  anal  vs  is .  also  known  as  the  "s+.-Hp  me+.Vwi"  n-r 
"cross  flow  hypothesis",  is  thus  an  approximate  one  in  the  sense  that 
a  certain  degree  of  interaction  between  adjacent  sections  is  neglected. 

The  cross  section  of  the  ship  at  x  will  now  be  taken  as  semi¬ 
circular;  the  correction  necessary  to  represent  other  ship  sections 
will  be  introduced  later.  Following  FJi.  Lewis  (15*29)  and  Weinblut  and 
St.  Denis  (1950),  the  flow  about  the  semi -submerged  body  used  in  the 
basic  derivation  will  be  assumed  to  be  identical  with  that  about  ">he 
lower  half  of  a  fully  submerged  body.  Corrections  to  account  for  the 
presence  of  the  free  water  surface  will  be  brought  in  later. 
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In  considering  the  pitching  and  hearing  notions  of  the  body  it  is 
necessary  to  introduce  a  second  coordinate  system  moving  with  the  ship 
with  its  origin  at  the  center  cf  gravity  of  the  ship.  The  longitudinal 
distance  of  any  section  of  the  ship  fro*  the  origin  is  designated  c 
(positive  forward).  Vertical  displacement  of  the  C.G.  (i.e.  the  heave) 
is  designated  by  s  (positive  upwards)  and  angular  displacement  or 
pitching  motion  is  designated  by  0  (positive  for  bow  displaced  upwards). 
The  vertical  displacement  of  the  section  at  x  due  to  pitching  is  then 
J;  Q  t  with  0  in  radians,  for  the  relatively  snail  angles  encountered. 
(It  is  also  assumed  that  cos  0^1.) 

The  two-dimensional  flow  pattern  between  the  control  planes  at  x 

1.  Vertical  velocity  of  the  center  of  the  circle 

T«t+|9-V0.  (6) 

2.  Vertical  component  of  wave  orbital  velocity 

•  2*yA  -2  k  h  c  2*yA  2*  ,  ^  v 

Tw  e  “  - - -  e  31  cos  ^  (x  -  ct)  (9) 

where  h  is  wave  amplitude,  X  is  wave  length,  v  -  -R  cos d. , 
the  depth  below  the  still -water  level  to  any  point  in  the  fluid, 
and  V7  ”  h  sin  2n  (x  -  ct)A» 

t 

3.  Apparent  variation  of  the  radius  r  of  the  ship  section 
at  the  control  planes  with  time;  r  *  r(t). 


All  motions  are  assumed  to  be  sufficiently  small  so  that  the  deri¬ 
vatives  of  the  potential  can  be  taken  on  the  surface  of  thj  circle  as 
if  its  center  were  at  its  initial  position  y  -  0,  and  the  known  expres¬ 
sion  for  the  potential  in  a  uniform  fluid  stream  can  be  applied,  despite 
the  slight  non-uniformity  induced  by  the  waves  which  are  assumed  to  be 
small. 


The  vertical  hydrodynamic  force  acting  on  the  length  dx  of  the 


submerged  semi-cylinder  is  given  by 

9/2 


f  =  2r 

dx 


p  cosoC  dd. 


(1< 


where  F  denotes  the  vertical  force,  is  the  polar  angle  as  defined 


in  Fi®.  il..  t  ire  -depend™  t  pnrt  of  the  prooour® 


c  f  py  0.0  UJWOJJICU 


from  Bernoulli's  equation.  Neglecting  the  squares  of  small  perturbation 


velocities 


p.pli 
p  \  at 


(i] 


where  0  is  the  velocity  potential  and  ^  the  mass  density. 


The  velocity  potential  of  the  flow  about  a  cylinder  due  to  the 


relative  vertical  velocity  (v  -v  }  is  given  by 


-(v  -  v  )  ^  cos  OC  . 


(12 


The  first  term  of  (12)  may  bo  considered  as  the  potential  due  to  the 
body  motion  in  smooth  water,  designated  0,  , 


The  second  tern  of  (12)  is  the  potential  due  to  body-wave  inter¬ 


action-  fit  . 

'  '  OW' 

„  r2 
*b»  ■  T»  -s  CMI*  • 

p 

2s  h  c  r  2syA  .  2s  /  -\  /.. 

■ - e  J/  coso(  cos  —  [i  -  ct).  (1 

The  potential  due  to  wave  notion  alone  is  0^, 

0 I  ■  h  c  $2«yA  cos  ^  (x  -  ct).  (1 

W  *v 

The  total  velocity  potential  is  the  sum  of  Equations  (13) »  (lii)  and  (15) 
*  "  *fce.  +  *bw  *  *v*  0 

Exciting  Forces 

Attention  will  now  be  concentrated  on  the  second  and  third  terras 
of  Equation  (16),  the  two  component  parts  of  the  velocity  potential 
which  give  the  exciting  forces  due  to  waves. 

The  pressure  due  to  is,  from  (ll)  and  (lh), 

_  k  s2c2h  n  r2  (-2s  R  cosa  )A  -  2n  A \ 

i>b.  ■  ?  ve - ^ —  °  "5  •  cosc'  8in  V  (l  * ct) 

-  ¥  -  ?  •(‘2’'  B  COi"*  )A  C»«  COB  (I  -  Ot). 

2 

On  the  surface  of  the  body  where  R  *  r  (and  since  c  "  g  X/2n), 
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f  2k 

X 


h  sin  ~  (x  -  ct)  r  cos<x.  e 


,(-2n  r  cos  of,  )A 


2  o  g  h  _  2it 


— cos  -T-  (x  -  ct)  £  coscf.  e 
'c  X 


(-2n  r  cos'*.)  A  (l?) 


(In  the  1955  paper  an  error  was  made  in  substituting  R  *  r  before  dif¬ 
ferentiating  with  respect  to  time,  since  r  is  a  function  of  time.) 


The  corresponding  component  of  the  force  due  to  body-wave  interfer¬ 
ence  is  obtained  by  substitution  of  (17)  in  Equation  (10 ): 

(f"|)bw  "  ^  g  h  r~  8ln  "T  ”  ct'  ”  ^  tLhJLZ.  cos  ^  U  -  ct) 


H 

lj*  COS20( 


(-2n  r  coscf  )/\ 


The  series  expansion  of  the  exponential  is 

^(-2n  r  coso?  ) A  ,  2nr  ,  2n2r2cos2o( 

a  *  j.  —  ~ t * ■  cos o\  ♦  ~  *  •  •  • 

K  x^ 


and  the  integral  is  evaluated  as 


„_2  (-2n  r  cost*  )A  ^ 

cos  o  d  oc 


n  _  k  n  r  +  3  ft^r2 
^  ^  8  X2 


Neglecting  cubes  and  higher  powers  of  the  small  quantity  rA  in  the 

2  2 

first  term  and  also  r  A  in  the  coefficient  of  the  t  term. 


/  p  2  2 

,d  ?n  „  ,  ,k  r  8  n  r  \  ,  2j  /  ^ 

(d^bw  ~2^hrj(2x- -J^r}  sin  t  (x  - ct) 

J  fr  /u  8n  r-i  2n  /  ,■>  ) 

+  c  (2  “  “ 5  C°S  “X  (x  -  ct)  f  * 

The  pressure  due  to  the  potential  of  wave  motion  0v  is. 
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from  (11)  and  (15), 

pv  •  ?  #  -  ^4^  -2”7A  -  T  (*- «*)■ 

At  the  surface  of  the  body, 

^  -  0  g  h  e(-2”  r  C08W  )A  sin  (x  -  ct).  (15) 

The  corresponding  component  of  the  force  is  obtained  by  substituting 


(19)  in  (10) t 


(f-|)w“  2  ^  g  h  r  sin  ^  (x  -  ct) 


n 

'2 


.  (-2n  r  coso()/l., 

coboI.  e  '  dot. 


where  the  integral  is  equal  to 


2  .22 

,  nr  .  h  nr 

1  ti  — r  * 

*  3  X 


Again  neglecting  cubes  and  higher  powers  of  r/X  t 


-2^>ghr(l-  +  j  ::-y-)  sin  ^  (x  -  ct).  (20) 

The  first  tern  of  Equation  (20)  is  seen  to  be  the  displacement 
force  resulting  from  the  wave  rise  or  fail  and  the  accompanying  increase 
or  decrease  of  volume.  The  second  and  third  terms  represent  a  modifica¬ 
tion  of  this  force,  due  to  the  ( approximate )  exponential  variation  of 
pressure  with  depth,  and  this  modification  is  known  as  the  "Smith  Effect" • 
The  entire  equation  represents  the  force  acting  under  what  is  usually 


referred  to  as  the  "Froude-Kriloff  Hypothesis".  This  force  which  is 
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exerted  by  the  waves  on  the  body  is  reduced  by  the  body*<wave  interference 
eiiect  indicated  by  Equation  (l8)j  for  a  send -cylindrical  section  at 
zero  speed  this  amounts  to  approximately  doubling  the  "Smith  Effect". 

Designating  by  p  the  angle  between  the  longitudinal  tangent  to 
the  body  surface  and  the  positive  x-axis,  the  derivative  r  is  evaluated 

>  ■  3J  H  ■  -v  tan  p  •  (211 

(In  the  1955  paper  the  sign  was  erroneously  taken  as  positive.)  With 
the  substitution  of  (2l),  the  sum  of  Equations  (18)  and  (20),  which  is 
the  total  exciting  force,  bee ernes 


12  -  (12)  ♦  (12) 

d  x  a  x  w  d  x  hr 


2p  g 


.  h  r  (  Cl  -  £  *  sl„  H  (x  -  et) 


.  V  A  „  ,n  8  n  r\  2n  , 

+  -  tan  p  (^ - j-j-)  005  "X  "  ct' 


This  expression  replaces  Equation  26  of  the  1955  paper.  It  differs  from 
it  in  the  sign  of  the  velocity-dependent  terms,  which  are  small,  and  in 
the  value  of  the  coefficient  of  the  (r  tan  p)/\  term,  which  is  also 
small. 

Two  steps  remain  to  be  taken.  First,  Equation  (22)  must  be  general¬ 
ized  for  ship  sections  other  than  semi-circular  and  second,  corrections 


must  be  introduced  for  free  surface  effects. 
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It  is  clear  that  the  factor  2  Q  g  t  h  sin  2*  (x  -  ct)/k  *  2  ©  g  r  -yj 
represents  the  change  in  displacement  force  with  wave  rise  and  fall;  r 
in  this  case  then  is  to  be  taken  as  the  half -be  aa,  B/2,  at  the  load  water¬ 
line.  Next  it  is  noted  that  when  the  body-wave  interaction  is  taken  into 
account  the  "Smith  Effect"  on  a  circular  body  is  doubled.  This  factor  of 
2  may  be  interpreted  as  (1  *  1c,)  by  analogy  with  G.l. Taylor's  expression 
(1928)  for  the  force  acting  on  a  body  placed  in  a  fluid  flow  with  a  velo¬ 
city  gradient.  Here  is  the  coefficient  of  accession  to  inertia  in 
vertical  flow  and  is  equal  to  1  for  a  circular  section. 

With  a  free  water  surface  and  the  formation  of  a  standing  wave 
system,  the  value  of  kg  “  1  for  the  circular  cylinder  is  modified  by 
a  factor  Milieu  is  designated  as  K^.  Ursell  (195U,  and  answer  to  Discus¬ 
sion,  Ko rvin-X.ro ukovsky ,  1955b)  has  computed  the  following  values  of  k^ 

2  2 

versus  urr/R  (or  u)  B/?g)  for  ths  circular  eylindsri  1 


u)2  B/2g  \ 


0 

©o 

.262 

.818 

.52k 

.632 

.785 

.592 

1.571 

.673 

2.09b 

.738 

2.356 

.762 

3.1U2 

.818 

3.927 

.85  9 

U.712 

.883 

In  the  absence  of  more  complete  information  it  will  be  assumed  that 


this  table  of  corrections  applies  to  non— circular  sections  as  well 
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Fro  a  experiments  with  an  oscillator,  Golovato  (1956).  derived  the 
coefficients  of  added  (virtual)  mass  in  heaving  oscillation  for  a  ship 
fora  symmetrical  fore-and-aft  with  U-sections  almost  wall-sided  at  the 
load  waterline.  Fig.  8  of  that  reference  shows  a  curve  very  similar  ■‘In 
trend  to  the  coefficients  of  the  above  table  but  with  values  about  20 % 
higher  and  with  the  minima  shifted  to  a  somewhat  higher  frequency.  The 
effect  of  these  differences  on  the  ship  response  are  expected  to  be  small. 

The  factor  (1  ♦  k^)c^)/2  will  then  be  applied  to  all  terms  of 
Equation  (22)  after  the  first  displacement  force  term.  In  the  earlier 
paper,  was  estimated  for  the  ship  as  a  whole  and  it  was  applied 

only  to  the  integrated  virtual  mass  and  inertia  effects  due  he  body 
motion  in  *wv>+.h  w  *+••»-?- •  ji+.  v«  crittcd  in  the  expression.  fc ~  the  exci¬ 
ting  forces.  Subsequently,  it  was  found  necessary  to  apply  the  k, 
correction  to  each  section  for  the  calculation  of  bending  moments  and 
highly  advisable  to  adopt  this  more  accurate  procedure  for  the  motion 
calculations.  This  omission  has  therefore  been  corrected  in  the  pre¬ 
sent  paper. 

Since  the  modified  "Soith  Effect"  terms  are  connected  with  virtual 
masses  and  since  the  effect  of  section  sh«pe  is  defined  by  the  coeffi¬ 
cient  k the  factor  r  in  this  case  is  interpreted  as  a  measure  of 


sectional  areai  i.  e 


r  * 


1/1 2  S/x 


(23) 
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and 


s 


where  the  limits  of  integration  s  and  b  are  the  values  of  x  at 
the  stern  and  the  bow,  respectively.  The  second  term  of  Equation  (27) 
results  from  the  consideration  that  the  water  pressure  acts  normal  to 
the  body  surface}  in  the  case  of  a  body  of  varying  circular  section 
the  moment  arm  is  (  §  ♦  r  tan  p)  and,  by  the  -use  of  the  relationships 
(23)  and  (2li),  for  a  norms!  ship  fora  the  moment  ana  may  be  assumed  to 
be  (§  +  d  S/n  d§  ).  Equations  (25)#  (26)  and  (27)  replace  Equations 
(33)  and  ( 3U)  of  the  earlier  paper. 

The  integrals  of  (26)  and  (27)  can  be  evaluated  readily  by 
Simpson's  rule.  By  changing  the  ship's  position  relative  to  the  wave 
the  maximum  values  or  amplitudes  of  the  exciting  force  and  moment  can 
be  found  as  well  as  the  phase  lags  C  and  T  .  Calculations  of  these 
amplitudes,  Fq  and  Mq,  were  made  for  a  5-ft-long  model  of  the  Series 
60,  0*60  block  coefficient  hull  (ETT  Model  No.  Uiii5)  in  waves  of  ship 
length  (\/L  ■  l)  and  wave  height  of  1.5  in.,  for  comparison  with  the 
experimental  data  described  in  lppendix  2  of  the  1955  paper.  A  very 
good  agreement  between  calculated  and  experimental  values  was  obtained 
except  at  zero  and  very  low  forward  speeds  (see  Fig.  15).  It  has  been 
observed  in  recent  years  that  at  such  speeds  there  is  often  interference 
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from  the  waves  reflected  by  the  walls  of  the  towing  tank,  while  at 
greater  speeds  this  interference  is  no  longer  encountered  and  the 
experimental  data  becone  more  reliable. 

Forces  Due  to  Body  Motions 

The  pressure  in  the  fluid  due  to  the  body's  own  notion  is,  from 
Equations  (ll)  and  (13) , 

Pbm  '  ?  TT2  "  -  R  *  h  c0SC*  '  C0S<X‘ 

■  -  v  r  co8C<.  ■*  2  ^  v  r  coso(  (28) 

at  the  body  surface  where  R  -  r.  The  distribution  of  vertical  forces  is 
obtained  by  substituting  (28)  in  (lO)i 

( jf— } .  =  —  f~i  ~  v  —  O  II  i>  ?•  tr  .  t  OO  \ 

<i  X'  DU  n  £  \  ‘ 

The  apparent  vertical  velocity  v  of  the  center  of  the  circular 
section,  given  by  Equation  (8),  consists  of  three  parts 
v  -  Z+J9-VQ, 

where  the  first  is  the  heaving  velocity,  the  second  is  the  vertical 

velocity  due  to  the  angular  velocity  of  pitching  and  the  third  is  the 

vertical  velocity  due  to  the  instantaneous  angle  of  trim  0  at  the 

control  planea.  Since  §  is  a  function  of  tires  and  §  *  -V, 

•  •  •  •  •  • 

v  -  •?.  +  ^  S  -  2  i  8  .  (30) 

After  substituting  (8),  (21),  and  (30)  in  Equation  (29),  the  following 


set  of  six  terns  is  obtained t 


Equation  (3l)  replaces  Equation  (39)  of  the  1955  paper.  Terns  (l) 
and  (2)  are  identical  with  (3)  and  (6)  of  Equations  (39)  and  tern  (3)  is 
the  sees  of  the  earlier  (l)  and  (5)  with  the  sign  of  the  latter  corrected, 
items  (k) i  (5)  and  (6)  are  twice  teres  (k)»  (7)  and  (2)  respectively  of 
the  earliar  Equation  (39),  which  were  incorrect  because  substitution  of 
R  ■  r  had  been  made  before  differentiation  with  respect  to  time. 


■j 

The  factor  (^  a  r  /2)  in  the  first  three  terms  of  (31)  is  evidently 

the  virtual  mass  of  an  element  of  body  length,  and  by  introducing  and 

k^  on  the  basis  of  the  reasoning  outlined  in  connection  with  exciting 

forces,  it  can  be  expressed  as  (  ^  S  The  factor  (^w  r  tan  £3) 

in  the  last  three  terms  is  the  derivative  with  respect  to  5  of 
2 

(^  n  r  /2)  and  so  it  can  be  expressed  as  d(^>  S  kgk^)/  dg,  . 


The  total 


force  due  to  the  body's  own  notion  is  then 


?ta  *  -S  /  S  Vli  (**£S-27  9)  d| 


/d  S  k-k, 

*  VPJ  - (*  *  5  0  -  V  0)  dg 


and  the  ncment  is 


“  "f  /s  Vt  (i‘ +  §  ®  - 2  7  *>§  d§ 


/d  S  k-k, 

#v?/“ds  (i*  §«-V0)s  d£ 


(32) 


where  the  integration  is  carried  over  the  length  of  the  hull. 


Since  and  are  functions  of  z  and  6  sad  their  deriva¬ 

tives,  the  terns  of  Equations  (32)  and  (33)  aay  be  transposed  to  the 
left-hand  side  of  the  coupled  equations  of  notion,  fbeir  contribu¬ 
tions  to  the  coefficients  of  these  equations  are  designated  by  the 


subscript  1.  Thus 
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/d(S 
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For  the  case  of  a  half-immorsed  spheroid  under  the  condition  of  a 
free  water  surface  but  neglecting  wave-making,  Havelock  gives  the  dyna¬ 
mic  coupling  terms  as 

-  p  M  V  0  for  the  heaving  force 
and  *  q  M  V  l  for  the  pitching  moment, 

where  M  is  the  mass  of  displaced  water.  In  general  p  f  q,  and  each 
is  given  by  a  fairly  complicated  expression  in  terms  of  ellipsoidal 
coordinates  and  associated  Legendre  functions  of  the  second  kind.  For 
a  long  SDheroid  Havelock  finds  that  p  ■  q  -  (l  ♦  k.  )/2,  or  .515  for  a 

X 

length-diameter  ratio  of  8  which  is  the  fineness  ratio  of  the  usual  ship 
form.  (Has kind  (19U6)  had  found  that  p  ■  q  for  a  thin  or  "Michell" 
ship  symmetrical  fore-and-aft.) 
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Thus  p  »  q  -  .  For  a  circular  section  =  1  .  and,  at  the 

oscillating  frequency  in  the  vicinity  of  synchronism  for  aost  ship  forms, 
is  of  the  order  of  0« V5«  It  appears  from  the  application  to  the 
models  in  this  paper  that  tes  damping  in  heave  is  reduced  ana  in  pitch 
is  increased  by  the  addition  of  the  dynaaic  coupling  terns. 

Dissipative  Daaping 

It  was  mentioned  previously  that  the  free  water  surface  was  not 
taken  into  account  in  the  basic  derivation  of  the  present  paper  and  that 
a  correction  for  it  must  he  introduced  independently.  The  effect  of  the 
free  surface  on  the  virtual  mass  has  been  allowed  for  approximately  by 
the  use  of  the  coefficient  k^  derived  by  Ursell  for  a  semi-cylinder 
(195U,  and  Discussion,  Xorvin-Kroukovsky,  1955b).  (Grim,  1953*  has 
also  calculated  this  effect  for  same  ship-like  forms  but  his  material 
is  not  extensive  enough  for  general  application.)  The  other  well  known 
effect  of  the  free  surface  is  the  dissipation  of  energy  in  the  formation 
of  waves  which  propagate  away  from  the  ship  in  all  directions.  In  the 
"strip*  method  of  analysis  it  is  assumed  that  waves  from  each  length 
segment  d  i;  propagate  laterally.  If  the  ratio  of  the  ®aplitude  of 
these  waves  to  the  amplitude  of  the  heaving  motion  of  a  ship  section  is 
designated  by  A  ,  the  damping  force  per  unit  vertical  velocity  v  of 
the  shipj segment  is  expressed  as  (Holstein  1936,  1937a,  1937b,  and 
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where  a)  Is  the  frequency  of  the  wares  radiated  '07  the  ship  and  is 
equal  to  the  frequency  of  ware  encounter. 


Holstein  (1936,  1937a  and  I937o)  and 


(19 k2'j  represented 


the  heaving  body  by  a  distribution  of  harmonically  pulsating  sources 
along  the  bottom.  In  that  case  the  asplitude  ratio  A  is  given  by 

A  *  2  e  *0^  sin  (k^y)  ( 

where  k  *  it)  c/g,  y  is  the  half -be  an  B/2,  and  f  *  S/5,  the  zxsz. 

O  V 

draft  of  a  ship  section.  This  theoretical  result  was  eouflmed  by 

Hcl  S  tC  i-H  1  3  (Teu  CApuTiM&UWC  j  ^rij£  X^3~wJ."tS  C  ^  wJC  »  T-  TtS^T;  '^4?  »P 

consistent,  yet  some  doubt  nay  be  felt  because  of  the  smallness  of  toe 
test  tank.  The  theory  is  approximate  and  acceptance  of  it  necessarily 
hinges  on  agreement  with  experimental  results. 


A  more  nearly  exact  theory  was  developed  by  Irsell  (l®h9,  1953, 
195k)  for  a  bearing  seaicylincer,  and  by  Iris  (195-3)  far  a  number  of 
analytically  defined  sections  closely  approaching  practical  ship  sec¬ 
tions.  In  the  case  of  a  sesi circular  section  Grim's  results  agree  with 
“-rseii5s.  In  their  tneory  the  damping,  force  is  calculated  as  a  boundary 


value  rrchlen. 


:rtunately  nc  emeriaentsl  verification  has  been 
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The  computed  coefficients  are  based  on  the  "strip"  method  of 


analysis  which  assumes  two-dimensional  fluid  flow,  while  Golovato's  and 
Haskind's  and  Riman's  tests  were  made  with  ship  models,  i.e.  in  three- 
dimensional  flow.  The  relation  of  damping  in  three-diiransional  flow  to 
dsnping  in  two-dimensional  flow  can  be  estimated  on  the  basis  of  the 
work  of  Havelock  (1956)  and  of  Vossers  (1956).  Havelock  calculated  the 
damping  coefficient  by  two  methods,  three-dimensional  and  "strip",  for 
a  submerged  spheroid  with  length-beam  ratio  of  8,  a  fair  value  for  com¬ 
parison  with  ship  models.  Vossers  made  similar  computations  from 
Has kind's  theoretical  results  (19U6)  for  a  "thin  ship"  in  the  sense 
of  the  approximations  introduced  by  Michell  (1896)  in  his  theory  cf 
wavs  resistance  of  ships.  The  results  of  both  are  reproduced  here  in 
Figs.  16,  17  as  ratios  of  the  coefficients  of  damping  in  heave  or 
pitch  by  the  two  methods. ^  It  should  be  remembered  that  accurate 

evaluation  of  damping  is  most  important  in  the  vicinity  of  synchro- 
nism.  The  values  of  the  frequency  parameter  u)  L/g  ( (y  L/g  in  the 
figures)  at  synchronism  are  given  in  the  last  two  columns  of  Table  5 
for  the  eight  models  to  which  the  computational  method  outlined  here 
has  been  applied.  (Also  included  in  Table  5  are  the  parameters  L/B, 
kg  and  k'  for  comparison  with  Havelock's  submerged  spheroid.) 


10 


The  subscript  S  in  the  figures  denotes  "strip"  method. 
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Omitting  from  consideration  the  original  trawler  and  yacht,  the  lowest 
values  of  the  parameter  u)  L/g  are  10.U  in  heave  for  the  T-2  tanker 
and  11.7  in  pitch  for  the  Series  60  hull,  while  for  the  models  with  dis- 
placement-length  ratio  of  60  ui  L/g  is  above  17*6.  At  these  valuse 
Figs.  16  and  17  indicate  negligible  corrections  for  damping  in  heave 
and  small  corrections  for  damping  in  pitch  (0  to  15£  based  on  Havelock’s 
computations  and  10  to  205t  from  Vossers’).  Since  the  corrections  for 
three-dimensional  effect  are  negligible  or  small  and  since  they  have  not 
yet  been  developed  for  ships  of  normal  foimt,  they  will  be  ignored  in  the 
present  work;  the  damping  will  be  taken  as;  computed  for  two-dimensional 
flow  by  Grim.  It  is  gratifying  nevertheless  to  have  measures  of  possible 
error  such  as  Figs.  16  and  1?  to  replace  arbitrary  assumptions  made  in 
the  earlier  paper  as  to  the  applicability  of  the  "strip”  method  of 
analysis. 

the  damping  force  for  each  section  is 
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The  contributions  of  these  expressions  to  the  coefficients  of  the 
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force  and  moment  equations  respectively,  depend  only  on  the  changes  in 
the  displacement  of  the  ship.  With  the  linearizing  assumption  these  are 
evaluated  on  the  basis  of  the  beam  of  a  section,  g  and  G,  the  coeffi¬ 
cients  of  the  angular  displacement  0  ,  depend  on  displacement  changes 
and  also  on  the  kinematics  of  the  fluid  flow  resulting  from  the  ship 
being  at  an  instantaneous  angle  of  trim. 

It  is  seen  that  the  damping  force  coefficient  b  in  heave  is  a 

function  of  frequency  of  encounter  but  is  independent  of  forward 

© 

speed  V  per  se.  This  appears  to  be  confirmed  by  the  experimental 
work  of  Gclovato  (1956).  However,  the  damping  moment  coefficient  B 
in  pitch  and  the  cross-coupling  -efficients  e  and  E  are  composed  cf 
dynamic  terms  proportional  to  V  and  dissipative  terms  independent  of 
V  (except  as  the  frequency  of  encounter  c<)  is  a  function  of  V). 

While  the  dynamic  terms  contribute  only  a  little  to  B,  they  make  very 
important  contributions  to  the  cross-coupling  coefficients. 
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TABLE  1 

MODEL  PROPERTIES 


Model  Numbers  ' 


lbb5 

1616 

litbb 

1723 

1699A 

1699B 

1699C 

1699D 

Length  BP,  ft. 

5.00 

5.00 

b.79 

5.71 

b.17 

- 

5.55 

- 

Load  WL,  ft. 

- 

- 

- 

- 

b.ll 

3.9b 

- 

- 

Rated  WL*,  ft. 

- 

- 

- 

- 

U.28 

U.28 

5.71 

5.71 

Beam,  ft. 

0.667 

0.667 

0.65 

0.608 

0.7b 

1.22 

0.555 

0.709 

Draft,  ft. 

0.267 

0.267 

0.29c 

0.208 

0.330 

0.690 

0.255 

0.U17 

Displacement(FW),  lb. 

J 

33.3 

33.3 

Ui.o 

2b.  5 

33.3 

50.8 

2b.  5 

2b.  5 

Block  Coefficient 

0.60 

0,60 

0.7b 

0.55 

0.51 

0,23 

0.51 

0,23 

A/(.01L)3 

122 

122 

171 

60 

20b 

370 

60 

60 

Radius  of  Gyration 
in  air,  ft. 

1.27 

1.27 

1.15 

1.37 

1.07 

1.07 

1.37 

1.37 

Natural  Pitching  Period 
afloat  in  calm  water, 

/ 

sec.  0.72 

0.65 

0.70 

0.59 

0.7b 

0.81 

0.59 

0.63 

Natural  Heaving  Period, 
afloat  in  calm  water. 

sec.  0.70 

0.70 

0.75 

0.6l 

0.73 

0.78 

0.61 

0.60 

* 

By  Cruising  Club  Rulei 

Rated  LWL  ■ 

0.3  WL 

♦  0.7 

,  where  WL  is 

the  load 

WL 

length  and  WL^0^  is  the  length  on  the  WL  at  a  draft  of  l.Ob  x  load  draft. 
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TABLE  2 

VALUES  0?  COEFFICIENTS  AND  EXCITING  F02CES  0? 
EQUATIONS  OF  MOTION  IN  HEAVING  ABC  PITCHING 
FOB  ETT  MODEL  HO.  11*1*5  in  WAVS  LENGTH  =*  SHIP  LENGTH 


0 

-17 


c 

0 

C.+  c. 

- 

0 

2 

x  2 

C 

136 

138 

d  *  D 

-.05 

-.05 

e2  '  E2 

-2.3 

-2.2 

*1  “E1 

0 

-.9 

e 

11 

CM 

1 

-3.1 

E 

C7 

-2.3 

- 

lii  7 
136 


X 

137 

-.05 

-2.1 

-1.6 

-3.7 

-jv  s 


-6 

-23 


-o 

-26 

-17 


-1 

185 

-.05 

-2.0 

-2.!i 

-Hill 

O.ii 


-10 

-27 


Values  of  parameters  obtained 

at  model  speed  (in  fpi 

Symbols 

0 

1.25 

2.25 

3.25 

iuOO 

£j 

6. 36 

7.93 

9.19 

10.I15 

11.2*2 

a 

1.37 

1.76 

1.72* 

1  H  (*. 
x*  »  '•> 

1.77 

A 

2.U0 

2.28 

2.23 

2.22 

2.22 

b 

5.6 

It.  7 

3.9 

3.0 

2.3 

B 

7.7 

6.3 

6.2 

5.2 

lut 

-V 
181 
—.06 
-1.3 
-3.0 
— *3 
1  0 


-9 

-26 


«F 

1.0 

1.0 

1.0 

0.3 

0.3 

*5 

3.6i 

3.7i 

3.7i 

3.6i 

3.51 

Based  on  uJt  *  0  when  wave  crest  is  at  station  11.  The  correction  of 
13°  is  to  be  applied  to  phase  lags  to  convert  to  aJt  *  C  at  station  10. 


BEST  AVAILABLE  COPY 
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TAELS  3 

7AI12S  3?  ZC&ZTZCI3FZS  SSL  2ICI73G  7CECZS  CF 
s^noss  of  xotios  3  $at3s  asi  pitchsc* 

POP.  377  MOI2L  50.  1616  3  ik*Z  LZSC-TS  -  SHIP  L2SC-IE 


" al'jes  of 

paraaetcrs 

obtained 

at  aodel  spe«c 

(.tJt  fpS 

Synbols 

0 

1.25 

2.23 

lL*  lO 

-4-i 

5.36 

~.?2 

?.l? 

n  i  r 

Ly- 

3.62 

a 

1.91 

1.77 

1.75 

_  •  t 

1.75 

A 

A  ''a 

-*  JV 

A  *  -A 

a  -»  — 

^  _ 

A 

;; 

c.ii 

b 

7.1 

5.? 

w.9 

L,  " 

3.6 

3 

11.1 

10.6 

9.7 

*  /• 

:-.c 

i  *  i 

c 

_ _  _ 

•  nn 

c 

_  -77X7 _ 

o 

p  p 

"1  '2 

0 

-1 

-k 

-3 

-12 

V 

227 

226 

223 

215 

215 

d  *  I) 

.015 

.015 

.015 

.016 

.013 

e2  “  *7 

-.1 

0 

0 

0 

0 

V  -*1 

0 

-.? 

—1.  6 

-2.1 

-3.0 

e 

-  .1 

-1.6 

-2.6 

-3.0 

3 

-.1 

.? 

1.6 

2.6 

3.0 

4- 

ao 

h 

0 

— 7.u 

-11.1 

-12.3  -13.6 

& 

5.2 

-2.2 

_x  o 

-6.1 

-6.6 

-  - _ 

__ 

>•  ^ 

? 

0.3 

0.3 

0.7 

0.6 

0.5 

V 

Lull 

6.i± 

6.11 

ii.  Oi 

3.91 
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TABLE  h 

CCXPJ&1SCE  Q?  SSlH-KiUHIlf*  ZEST3CTZB  KOZSL  acmOKS 

15  TWI5G  TJBTK  I2BEGCUE  SID  SgJkS, 

! 

AS  0BTAI52D  BT  TSffi  DI7722EK?  !«ET33DS 


Beam,  Inches _  _ Pitch,  degrees 

0  speed  2.53  ft/sec.  0  3pfegd  2.53  ft/scc 

3y  direct  analysis  of  aodel 

■otions  in  irregular  wares  1.1  1.3  5.2  6.7 

By  calculation  fro*  ware  spectra 
and  experlsaen  tally  aeasured  res¬ 
ponses  to  regular  wares  1.1  1.2  U.3  6*1 

By  calculation  fro*  ware  spectra 
aid  Hialytlcally  eomrated  res¬ 
ponses  to  regular  wares  1.1  1. 3  5.0  5-3 


*Mesn  of  1/10  hipest  aaplitrodes 
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!i=U  5 

TiLCES  ul^L/z  ?CE  33SCS3MCS  SPEED,  L/B,  c.  and  *•  OF  7$  7J2Z30S 
SEE?  FOatS  7CE  GXPAEISOS  tfTTE  ECBL0C17S  5C3S27 2D  SF9S2CID 


Model 

1/5 

*■ 

c. 

id* 

xl 

£ _ 

c 

Pitch 

Dim 

Pitch 

oeacre 

Hatelock’s  Sub—rged 

Spheroid 

•  ?fc5 

.5i 

- 

- 

- 

- 

Series  60  (Model  1LL<3) 

.72 

.«2 

t.  7 

?.o 

H.7 

y 

7 -bast  (Model  1616} 

•'  O 

•  Tt. 

?* ' 

?.c- 

— i-*  “ 

12-5 

7-2  lantaer  'Model  IcLe} 

?Ji 

.63 

•  51 

?.o 

3.L 

lt-0 

10^ 

Destroyer  (Model  i?23) 

?.is 

.90 

• 

10.6 

10.3 

20.1 

15.5 

Trawler  (Model  l69?A) 

.66 

3.c 

3.6 

9.5 

sr.Q' 

Tacit.  (Model  169?B) 

> 

.70 

.do 

•  •  ' 

".1 

;  .r 

w- 

Lecgtiieaed  Trader 
(Model  1699C) 

10.3 

e> 

.Oh 

.53 

10.6 

10.3 

20-1 

15.3 

lengthened  Tacit 
(Model  169^3) 

3.0 

.72 

.la5 

10.0 

10.5 

17.6 

Ip.ii 

,  katc  X 

i_  for  t£*  e=±d  Tarns  >xits  free  surface  correcbicc  i*  j  *  w  1  —y 
~  a.  Maes  at 


Displaced  bj  a 


Added  Xanect  of  Inertia 
Moser^t  of  Inertia  oT  aarcer 
Displaced  tj  s.  Snip. 
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Figs.  6  and  h  Substitcue  -e,  *  “or 


-*i '  -h  ( 


corrected  for 


Trans  actions  of  SK1W£) , 


i  Observed  data  on  pnase  relation¬ 
ships  ojaitttid  for  Hoeel  1616, 

X/L  *  1.50,  will  be  shewn  in 
SSiMS  'Trans actions. 


Fig.  5  Modal  I6S9B  ('TochO 


FIG.  ft  Motions  of  5~Ft.  E.T.T.  Mortal  s  1  JJj5  and  1  61  in  Waves  of 
1.0  and  1,$  Model  Lengths  by  1.2 5  in.  High.* 


Model  Speed,  Ft/Sec 


FIG.  9  Motions  of  U.8**Ft.  E.T.T.  Model  11*1*2;  in  Waves  1*.8  ft.  Long 
by  1.2  in.  High.*  j 


Phase 


500 


300 


c 

?2 


Double  AmpStudes 
-Model  1723.  X/L»(25- 

m 

1 

■ 

it: 

n 

n 

_z 

J?° 

55 

£°o  2  $  & 


Phase  Logs^^ 
"Model  1723,  X/L-125- 

HKm 

k^T 

_ - « 

r  _ 

■  rrrcrt 

r0- 

1 

Model  Speed,  FI /Sec 


El 

a.  a 


1  Doubks  Amplitudes 
J-Model  l723,X/L*2.0Cr 

—  1 

S 

i 

8 

tn 

J~D  J 

_ *_iJ 

,  j 

r” 

_ 

5  2 

* 6  8 

Phase  Logs 
hModel  i723VL«2.(Xrl 


HeawS 

' 

— — 

r 1  ■  1  "1 

* 

_ 

Pircrr 

100 

G 


t — §  2 — % — 5 — 8 

Model  Speed,  Ft/Sec. 


1U  Mo-oians  of  p.7-ft.  E.I.i.  Hocci  1723  in  waves  of  l.u,  1.23, 
1.5j  and  2.0  Model  Lengths  by  1*U3  in.  High.* 


Circles  indicate  E.T.T.  experimental  data  (open  for  heave,  solid  for 
pitch)  and  curves  show  calculated  notions.  Zero  phase  lag  corresponds 
to  maximum  pitch  or  heave  with  wave  nodal  point  at  midsection. 


Phase  Lag- Degrees  Phase  Lag- Degrees 


t  l _ U±  J _ I^J _ I _ 1 _ =*1 _ 1 _ j  '  i - j - I 

(02  46  80  2468  10  I*  14 

FrtOfRey  Of  EMOffiW.  U« 


FIG,  13  Motion  Spectra  Predicted  by  Three  Different  Methods  for 
Destroyer  Model  1723» 

1)  By  direct  analysis  of  motions  in  irregular  waves  (solid  lines) 

1 

2 )  By  ealcu5  +-ion  from  wave  spectra  sad  experimentally  measured 
responses  regular  waves  (dash  lines), 

3)  By  calculati.  or  wave  spectra  and  analytically  computed 
responses  to  regular  waves  (♦)• 


r:.G,  lii  Sketch  Illustrating  Notation  Used  in  Appendix. 


Forward  Speed, V,  Ft/9*c4A9cHnrt  The  Warn 


FIG.  15  Comparison  of  Conpmted  and  Ejqperiaen  tally  Measured  Exciting 
Force  and  Moment  Amplitudes  for  E.T.T.  Model  LUo  in  Mares 
5  Ft.  by  1.5  in. 


best  available  copy 


FIG.  16  Ratios  of  Three-Dimensional  to  Two-Dimensional  Calculations 
of  Damping  Coefficients  for  Submerged  Spheroid  of  L/B  *  8. 
(Fig.  1  of  Harolock,  1956) 


FIG.  17  Ratios  of  Three -Dimensional  to  Two-Dimensional  Calculations 
of  Damping  Coefficients  for  Thin  "Michell"  Ship  of  B/H  *  2. 
(Fig.  2  of  Vossers,  1956) 


